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We consider photonic crystal fibres (PCFs) made from arbitrary base materials and 
introduce a short-wavelength approximation which allows for a mapping of the Maxwell's 
equations onto a dimensionless eigenvalue equations which has the form of the Schrodinger 
equation in quantum mechanics. The mapping allows for an entire analytical solution of 
the dispersion problem which is in qualitative agreement with plane-wave simulations of 
the Maxwell's equations for large-mode area PCFs. We offer a new angle on the foundation 
of the endlessly single-mode property and show that PCFs are endlessly single mode for 
a normalized air-hole diameter smaller than ~ 0.42, independently of the base material. 
Finally, we show how the group-velocity dispersion relates simply to the geometry of the 
photonic crystal cladding. 
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1 Introduction 

Photonic crystal fibres (PCF) are a special class of optical fibres where a remarkable 
dielectric topology facilitates unique optical properties. Guiding of light is accomplished 
by a highly regular array of air holes running along the full length of the fibre, see Fig. ^ 
In their simplest form the fibres employ a single dielectric base material (with dielectric 
function = ?t|) and for the majority of fabricated fibres silica has been the most common 
choice [1,2]. This preference has of course been highly motivated by silica's transparency 
in the visible to infrared regimes, but has also been strongly driven by the highly matured 
glass technology used in the fabrication of standard telecommunication fibres. However, 
a growing interest in fibre optics for other wavelengths and light sources has recently 
led to a renewed interest in other base materials including chalcogenide glass [3], lead 
silicate glass [4], telluride glass [5], bismuth glass [6], silver halide [7], teflon [8], and 
plastics/polymers [9]. 

The huge theoretical and numerical effort in the PCF community has been important in 
understanding the dispersion and modal properties, but naturally there has been a clear 
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emphasize on silica-based PCFs. PCFs fabricated from different base materials typically 
share the same overall geometry with the air holes arranged in a triangular lattice and 
the core defect being formed by the removal of a single air-hole. However, it is not yet 
fully clarified how the established understanding may be transferred to the development 
of PCFs made from other base materials. To put it simple; PCFs made from different 
base materials share the same topology, but to which extend do they also have optical 
characteristics in common? 

When addressing this question it is important to realize that the otherwise very useful 
concept of scale-invariance of Maxwell's equations [10] is of little direct use in this par- 
ticular context. Since PCFs made from different base materials do not relate to each 
other by a linear scaling of the dielectric function, e(r) s 2 e(r), scale invariance can- 
not be applied directly to generalize the results for e.g. silica to other base materials. 
Electro-magnetic perturbation theory is of course one possible route allowing us to take 
small changes of the base-material refractive index into account, but in this paper we will 
discuss and elaborate on an alternative approach which was put forward recently [11]. 

As a starting point we take the classical description based on Maxwell's equations, which 
for linear dielectric materials simplify to the following vectorial wave-equation for the 
electrical field [10] 

u 2 

V xV x E(r) = e(r)— E(r). (1) 

cr 

Here, e is the dielectric function which in our case takes the material values of either air or 
the base material. For a fibre geometry with rn along the fibre axis we have e(r) = e(rj_) 
and as usual we search for solutions of the plane- wave form e'^n ~ Uft ' with the goal of 
calculating the dispersion a>(/3) relation. 

2 Short-wavelength approximations 

Today, several approaches are available for solving Eq. (0) numerically, including plane- 
wave, multi-pole, and finite-element methods [12-14]. Such methods will typically be 
preferred in cases where highly quantitative correct results are called for. However, nu- 
merical studies do not always shine much light on the physics behind so here we will 
take advantage of several approximations which allow for analytical results in the short- 
wavelength limit A <C A. The key observation allowing for the approximations is that 
typically the base material has a dielectric function exceeding that of air significantly, 
Sf, 3> 1. In that case it is well-known that the short- wavelength regime is characterized by 
having the majority of the electrical field residing in the high- index base material while 
the fraction of electrical field in the air holes is vanishing. 

The, at first sight very crude, approximation is simply to neglect the small fraction of 
electrical power residing in the air-hole regions. Mathematically the approximation is 
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implemented by imposing the boundary condition that E is zero at the interfaces <9f2 to 
air. Since the displacement field D = eE is divergence free we have = e'V -E + E- Ve f=s 
e'V ■ E and the wave equation (JTJ) now reduces to 

0, (2) 



-V 2 E{r) ne b ^E(r), E(r) 
c z 



redn 

and since solutions are of the plane- wave form e*^ r n ~ ut > along the fiber axis we get 

> 2 

-;2 



-VlE(r ± )+f3 2 E(r ± )^e b ^E(r ± ), E(r ± ) = 0, (3) 



where V 2 ± is two-dimensional Laplacian in the transverse plane. The essentially scalar 
nature of this equation makes the problem look somewhat similar to the more traditional 
scalar approximation 

Vi* (r ± ) + e(r ± )^ ¥(rj « (3 2 9(r ± ), (4) 

that has been applied successfully to PCFs in the short-wavelength regime [15]. However, 
while that approach took the field and the dielectric function in the air holes into account 
we shall here solve Eq. © with the boundary condition that E is zero at the interfaces to 
air. We note that this of course is fully equivalent to solving Eq. (@J) with the boundary 
condition that \l/ is zero at the interfaces to air. Obviously, the scalar problem posed by 
Eq. (jnj) is separable and formally we have that 
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yjtt{ + Q 2 = ^Ql + (c(3/n b ) 2 (5) 



where Q» = cf3/n b is the frequency associated with the longitudinal plane- wave propaga- 
tion with a linear dispersion relation (i.e. fin oc (3) and Q± is the frequency associated with 
the transverse confinement/localization. At this point we already note how an arbitrary 
base-material refractive index n b enters in the frequency associated with the longitudinal 
plane-wave propagation. This is solely possible because we have a Dirichlet boundary 
condition in Eq. (JHJ) which does not depend on the refractive index of the base material. 



3 The Schrodinger equation eigenvalue problem 

In the following we rewrite Eq. (JSJ) as 

u= —J^A-z + pz (6) 
n b 

where now 7 is a dimensionless number characterizing the confinement/localization. It is 
of purely geometrical origin and thus only depends on the normalized air-hole diameter 
d/A. From Eq. (J3J) it follows that 7 is an eigenvalue governed by a scalar two-dimensional 
Schrodinger equation 

-A 2 V^(r ± )= 7 2 ^(r ± ), ij>{r ± ) =0. (7) 
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The same equation was recently studied in work on anti-dot lattices in a two-dimensional 
electron gas [16] and there will be many similarities between such quantum mechanical 
electron systems and the present electromagnetic problem. In the language of quantum 
mechanics the wavefunction if; is subject to hard-wall boundary conditions corresponding 
to an infinite potential barrier in the air-hole regions. The task of calculating the optical 
dispersion properties has now reduced to solving the scalar two-dimensional Schrodinger- 
like eigenvalue problem posed by Eq. (JJJ). The strength is obviously that Eq. (JJJ) is material 
independent and thus we may in principle solve it once and for all to get the geometry 
dependence of 7 and thereby the optical dispersion of the PCF. In the following we solve 
Eq. (jJJ) for different classes of modes and compare approximate results to corresponding 
numerically exact finite-element solutions. 



3.1 The fundamental space- filling mode 

The fundamental space-filling mode is a key concept in understanding the light guiding 
properties of PCFs. It is the fundamental de-localized mode of Eq. in an infinite 
periodic structure and the corresponding mode index corresponds to an effective material 
index of the artificial periodic photonic crystal lattice. In the language of band diagrams 
and Brillouin zones it is defined in the T-point where Bloch's theorem is particular simple; 
as illustrated in Fig. El the wave function is subject to a simple Neumann boundary 
condition on the edge dQ of the Wigner-Seitz cell, i.e. 

-A 2 V^(r ± ) = 7c >K), 4>(r ± ) = n ■ V^(r ± ) . = 0, (8) 



where n is a unit vector normal to the boundary dQ. We thus have to solve an eigenvalue 
problem on a doubly-connected domain. By a conformal mapping, leaving the Laplacian 
and the boundary conditions invariant, one may in principle transform the hexagonal 
shape with the circular hole of diameter d into a simple annular region of inner diameter d 
and outer diameter 2R [17] leaving us with an eigenvalue problem in the radial coordinate. 
Finding the exact conformal mapping may be a complicated task so here we simply use 
the approximate mapping with 




fl = At/-£-~0.53A (9) 



which serves to conserve the area of the Wigner-Seitz unit cell. Neglecting any small 
changes to the eigenvalue equation caused during the mapping we get the following eigen- 
value equation in the radial coordinate, 

-A 2 («9 2 + r- l d r )^{r) « 7c>(r), *i>(d/2) = ^P'(R) = 0. (10) 

The solution is given by a linear combination of the Bessel functions Jo and Yq and the 
eigenvalue is determined by the roots of the following expression 

^i^Hi^-^n^- < ii; 
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In general the equation cannot be solved exactly except for at the point 

as is easily verified by insertion and using that a n ^ m is the mth zero of the nth Bessel 
function, i.e. J n {o^n,m) = 0. In particular we have that a^i — 2.405 and a± t i ~ 3.8317. 
However, expanding the left and right-hand sides of Eq. (|12|) around this point to first 
order in both 7^ and d/A we get an equation from which we may isolate 7^ yielding an 
approximate analytical solution of the form 

M Ci + c^7a) 2 ' (13) 

The coefficients are given by expressions involving Bessel functions Yq and Y\ evaluated 
at «o,i and a^i, see Appendix, but for simplicity we here only give the corresponding 
numerical values of these expressions; C\ ~ —0.2326, C2 — 2.7040, and C3 ~ 1.0181. 



3.2 Core-defect modes 

While the fundamental space-filling mode, discussed above, is a de-localized mode it is 
also possible to have strongly localized states, especially in the presence of defects such 
as one formed by removing an air hole from the otherwise periodic lattice of air holes, see 
figures ^ and |2j In this way light can be guided along the defect thus forming the core of 
the fibre. The requirement for corralling the light is of course that the core-defect mode 
has an eigenvalue 7^ not exceeding the corresponding eigenvalue 7^ of the surrounding 
photonic crystal cladding. To put it slightly different the defect in the air-hole lattice 
effectively corresponds to a potential well of depth 7^ and a radius roughly given by 
p = A- d/2, see Fig. El 

In order to calculate the fundamental core-defect mode we first note that in the limit 
of d/A approaching 1, the problem can be approximated with that of a two-dimensional 
spherical infinite potential well with radius p = A — d/2, see Fig. |3J For this problem the 
lowest eigenvalue is 

7) = (TTIif (14) 

Although this expression yields the correct scaling with d/A, the approximation obviously 
breaks down for small values of d/A. In that limit we follow the ideas of Glazman et 
al. [18], who studied quantum conductance through narrow constrictions. The effective 
one-dimensional energy barrier for transmission through the constriction of width A — d 
between two neighbouring holes has a maximum value of 7r 2 in the limit where d goes to 
zero. We thus approximate the d/A — > problem with that of a two-dimensional spherical 
potential well of height 7r 2 and radius A, 

-A 2 (<9 2 + r~ l d r )ij(r) + 7r 2 9(r - A)^(r) = f Ctl tp(r), (15) 
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where 0(x) is the Heaviside step function which is zero for x < and unity for x > 0. 
The lowest eigenvalue 7^ for this problem is the first root of the following equation 



In the case where the potential 7r 2 is replaced by an infinite potential the first root is given 
by ojqp Lowering the confinement potential obviously shifts down the eigenvalue and in 
the present case the eigenvalue is roughly 7^ ~ ir. In fact, expanding the equation to 
first order in 7 Cj i around e.g. \fii we get a straight forward, but long, expression with the 
numerical value 7^ ~ 3.221 which is in excellent agreement with a numerical solution of 
the equation. Correcting for the low-rf/A behaviour we thus finally get 



7 c 2 i^7 c \+^l- lim \ a ° A 2 =7 c 2 i+ ^ A ^ agi- (17) 

For the first high-order mode we may apply a very similar analysis. This mode has a 
finite angular momentum of ±1 with a radial J\ solution yielding 

(a - A) A 

n, 2 ~^ 2 +AZ — A n 2 (IR) 

7c,2 — 7c,2 + , d ,2 "1,1' V iC V 

Here, 7 2 2 — 7.673 is the second eigenvalue of the two-dimensional spherical potential 
well of height ir 2 and radius A. Again, it can be found from an equation very similar to 
Eq. ffToT). 
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3.3 Two-dimensional finite-element solutions 

Above we have solved the geometrical eigenvalue problem analytically by means of var- 
ious approximations. In this section we compare the quality of these results by direct 
comparison to a numerically exact solution of the two-dimensional eigenvalue problem. 
The developments in computational physics and engineering have turned numerical solu- 
tions of partial differential equations in the direction of a standard task. Here, we employ 
a finite-element approach [19] to numerically solve Eq. (JJJ) and calculate 7 2 versus d/A. 
We employ an adaptive mesh algorithm to provide efficient convergence. For cladding 
mode we implement Eq. (jHJ directly while for the defect modes we solve Eq. ((2)) on a 
finite domain of approximate size 10A x 10A with the defect located in the center of 
the domain. For d/A down to around 0.1 we have found this to sufficient to adequate 
convergence for the strongly localized defect states, thus avoiding strong proximity from 
the domain-edge boundary condition which for simplicity has been of the Dirichlet type. 
Figure summarizes our numerical results for the first cladding mode 7^ as well for the 
two first core-defect modes 7 2 x and 7 2 2 . The dashed lines indicate the corresponding 
approximate expressions in Eqs. (fT3*|) . (fTTj) . and (fTHj) obtained by analytical means. As 
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seen the qualitative agreement is excellent, but in order to facilitate also quantitative 
applications of the results we also include the results of numerical least-square error fits 
below. The thin solid lines show the expressions 



^-3.666+ \ A ;.£ og,i, (19a) 



7c 2 2 ^ 8.691 + ^ ^-f a\ x , (19b) 

(2-1) 

7 S ^ f-2.82476 + 5 - 23695 , /A ^| 2 , (19c) 

,d V 1.17908 - d/AJ ' v 7 

which match the numerical data within a relative error of less than 2% around the most 
important cut-off region d/A ~ 0.42. 



4 Derived fiber optical parameters 

With Eq. (jUJ) at hand we have now provided a unified theory of the dispersion relation in 
the short- wavelength regime for PCFs with arbitrary base materials and Eq. (jUJ) illustrates 
how geometrical confinement modifies the linear free-space dispersion relation. 

In fibre optics it is common to express the dispersion properties in terms of the effective 
index n c s = cf3/u versus the free-space wavelength A = c2n/u. From Eq. (P) it follows 
straightforwardly that 

n cS = nftWl - - \ , ( T ) (20) 



47T 2 7t| \A, 

which obviously is in qualitative agreement with the accepted view that n e g increases 
monotonously with decreasing wavelength and approaches in the asymptotic short- 
wavelength limit as reported for e.g. silica-based PCFs [2]. Similarly, the group- velocity 
v g = du/d/3 becomes 



V »=n b f-^ b {V (21) 
while the phase velocity v p = u / (3 becomes 

v p = - 1 (22) 

Y 1 47r2nf \a) 

The group-velocity dispersion is most often quantified by the wave-guide dispersion pa- 
rameter D wg = d(l/v g )/d\ which becomes 



D 1 ^ 



cA 47r 2 n{, 



An 2 nl \ A 



-3/2 

(23) 
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where any possible material dispersion nt(A) has been neglected. We note that since 
D wg > the geometry of the air-hole lattice always causes a positive wave-guide dispersion 
parameter. Large-mode area PCFs belong to the regime with A C A were we predict the 
following general magnitude of the wave-guide dispersion of a large-mode area PCF 

1 l\ 

lim D wg = -r-r^-. (24) 
a<a B cA47r 2 n 6 



Finally, the recently suggested parameter Vpcf = A^//3^ — (3^ [20] can be shown to be 
a purely geometrically defined parameter in the large- mode area limit. From Eq. © it 
follows straightforwardly that 



hmV PC F = ^l 2 cl -ll v (25) 

This implies that the endlessly single-mode property [2, 20-22] is a wave phenomena 
independent of the base material refractive index of purely geometrical origin. Higher- 
order modes are only supported for d/A > 0.42, as seen in Fig. for which Vpcf ^ n [20]. 



5 Comparison to fully-vectorial plane-wave simulations 

In the previous sections substantial analytical progress was made, but we still remain to 
address the question to which extend the basic assumption behind Eq. Q holds except 
that it becomes exact for the fundamental modes as A/A approaches zero. In Fig. 13 we 
compare the analytical results for the effective index in Eq. (|2T)|) to results obtained by 
fully-vectorial plane-wave simulations [12] of Eq. (0). For the fundamental space-filling 
mode we have employed a basis of 2 6 x 2 6 plane waves while for the core-defect modes we 
have employed a super-cell configuration of size 10A x 10A in the case of d/A = 0.4. 

Panel (a) shows results for the fundamental space-filling mode for various values of the 
base refractive index. As clearly seen the theory agrees excellently in the short- wavelength 
limit while pronounced deviations occur as the wavelength increases and the field pen- 
etrates deeper into the air hole regions and also vectorial effects become important as 
reported in [15]. Similar conclusion applies to the corresponding results in panel (c) for 
the fundamental core-defect mode. While d/A = 0.4 is a case of particular technological 
interest because of its endlessly single-mode property we would like to emphasize that 
equivalent agreement is found for other values of d/A (results not shown). However, since 
the success of the approximation in Eq. (jHJ) really is that A/n^ <C 2p = 2A — d and 
A/n^ <C A — d for the core and cladding results, respectively, there will of course be a 
small (^-dependence. However, this also indicates that the agreement will increase with 
an increasing refractive index rib of the base material as has also been observed for even 
higher values of rib than those studied in Fig. |H1 (results not shown). 
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6 Discussion and conclusion 



It is today more than 10 years ago that the PCF was invented by Russell and co-workers [1] 
and there is a still growing community of researchers directing their efforts toward fabri- 
cation and experimental studies of silica-based PCFs as well as quantitative modelling of 
the optical properties. At the same time the community is obviously facing the challenges 
and opportunities of new exciting fiber materials. PCFs made from different base mate- 
rials share the same topology so it seems quite natural to assume that they also have at 
least some basic optical properties in common. Let's return to the question posed in the 
introduction: to which extend do PCFs made from different materials have optical char- 
acteristics in common? The present work do to some extend address this question. Most 
importantly we illustrate how the waveguide dispersion originates from the geometrical 
transverse confinement /localization of the mode and how the endlessly single- mode prop- 
erty arises as a sole consequence of the geometry which acts as a modal sieve. In particular 
we have shown how PCFs are endlessly single-mode for d/A < 0.42 irrespectively of the 
base material an analytical expression for the wave-guide dispersion parameter applicable 
to large-mode fibres. 

For small-core PCFs our theory provides qualitative correct results though more quantita- 
tive insight still calls for fully vectorial simulations, see e.g. [23]. However, for large- mode 
area PCFs our expressions not only gives qualitative insight, but also quantitative correct 
expressions which may be used in straightforward design of large-mode area PCFs with 
special properties with respect to the group-velocity dispersion and the susceptibility of 
the fundamental mode with respect to longitudinal non-uniformities [24]. 
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A Taylor expansion of Eq. ( 1111) 

Taylor expanding the left and right-hand sides of Eq. (jllj) around the point 




(26) 
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to first order in both 7^ and d/ Awe get 

t 'd-d \ t Y?(a 1)X ) - y o 2 (ao,i) (id - 7o 



do J ^l 2 («i,i) V 7o 

^i(Q!o,i) d f d - d \ , , 

(7ci - 7o) 



Y (a ,i) A \ d 
+ O([d-d } 2 )+O([ lcl - l0 ] 2 )- (27) 

Next, solving for 7 c i we may write the result in the form of Eq. (jlHj) with 

r 1q(«o,i) ( ^oV 1 (0Q \ 

Ci=1o- 777 r -j- , (28a) 



Fi(a ,i) VA 

C2_ 7on 2 KOn 2 Ki) UJ ' (28b) 

r _ ^o 3 («o,i) ~ *o(<*o,i) >f («i,i) do r9S x 

3 " 7o^iKi)n 2 Ki) A' 1 j 

which have the numerical values listed in the paragraph below Eq. (|13|). 
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Figure 1: Optical micrograph of a silica-based PCF with a pitch A of order 10 microns 
and d/A ~ 0.5. Courtesy of Crystal Fibre A/S, www.crystal-fibre.com. 
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Figure 2: The topology of a PCF formed by omitting an air hole of diameter d in a 
triangular lattice of air holes with pitch A. Dirichlet boundary conditions are applied 
to the boundary dVL indicated by red solid lines. The red dashed line indicates twice an 
effective core radius p — A — d/2. 
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Figure 3: The hexagonal Wigner-Seitz unit cell of the periodic triangular lattice with 
lattice constant/pitch A. The dashed line indicates the annular region which has the same 
the doubly-connected domain of hexagonal outer shape. Dirichlet and Neumann 
boundary conditions apply to the inner and outer boundaries dQ and dfl, respectively. 




Figure 4: Zoom in on the core-defect in the extreme limit of d/A — 1. A Dirichlet bound- 
ary condition applies to the boundary dQ and the dashed circle indicates the approximate 
core radius p = A — d/2. 
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Figure 5: Geometrical eigenvalues 7 2 versus normalized air-hole diameter d/A. The data 
points are results of finite-element simulations while the corresponding dashed lines are 
the approximate results in Eqs. (|1HJ). (|17jl. and fjTHf) . The thin solid lines shows the 
numerical fits listed in Eqs. (j!9aj) . ()19b|) . and ()19c|) . Note how only a single core-defect 
mode is supported in the shaded region, d/A < 0.42. 
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Figure 6: Effective index versus normalized wavelength A/A for the fundamental space- 
filling mode (a) and the fundamental core-defect mode (c) for a PCF with d/A = 0.4 
with a base material refractive index n& varying from 1.45 to 1.6 in steps of 0.05 from 
below. The data points are the results of a vectorial plane-wave simulation of Eq. (0) while 
the dashed lines show the corresponding results based on Eq. (|2Uj) with Eqs. ()19a|) and 
(|19c|) for geometrical eigenvalues 7^ and 7^. Panels (b) and (d) show the corresponding 
eigenf unctions ifj. 
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